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i rings as an option for hydrogen targets. The change of the hydrogen content due to 

^ I the radiation dose applied by the circulating proton beam has been investigated in 
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that the loss of hydrogen as a function of applied dose receives contributions from 
a fast and a slow component. 
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1 Introduction 



Polymers are widely used as target materials in nuclear- and particle-physics 
experiment, since they provide an effective proton target due to the large hy- 
drogen content. As an example polypropylene fiber targets of 4 x 5 yum^ cross 
section have been used for the EDDA-expcrimcnt [1] at COSY [2] to mea- 
sure precise excitation functions of the unpolarized differential cross section 
at beam energies from 0.5 to 2.5 GeV. 

For polypropylene the H:C ratio should be 2:1, however, upon irradiation of 
the target this ratio is going to decrease due to so-called cross-linking [3,4] 
of the CH2-chains, a process in which H2-molecules are released when C-C 
bonds replace two C-H bonds. Losses of hydrogen or carbon due to nuclear 
reactions of the beam protons are negligible in comparison. 

In an internal target experiment both the beam position and width change 
during acceleration [5,6], such that an experiment will sample different areas 
of the target at different beam-momenta [6] . As the applied dose will also vary 
with beam position the effective thickness of the hydrogen target will become 
a function of beam momentum, as soon as the area of the beam-target overlap 
changes. 

In this paper we investigate the dependence of the hydrogen loss on the applied 
dose by irradiating a CH2 fiber target under very controlled conditions. The 
number density, i.e. hydrogen atoms per volume, for polypropylene is given by 
(see Tab. 1 for our notation) 

riHfi = — TT ■ (1) 



which will change when the target has been irradiated with some dose D. If 
we introduce a normalized hydrogen density ^Pat, it will be a function of the 
applied dose, and thus of position x along the fiber and time: 

gj,{x, t) = g^{D{x, t)) = I^^Ri^. (2) 



The aim of this paper is to determine qn{D). We used a stored proton beam at 
the COSY-accelerator at Jiilich on a CH2-fiber target, and recorded the ratio 
of the pp-elastic scattering yield to that of proton-carbon inelastic scattering 
as a function of time. This ratio directly determines dose related changes in g]\f. 
The scattering yield of elastic pp-scattering allows a precise determination of 
the luminosity and the COSY beam parameters, like width and position. Using 
this information the applied dose D{x,t) can be calculated with accuracy. 
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dy, dz 

A = dy ■ dz = 20 ^im^ 
Na 

Qcn, = 0.905^ 
^CH. = 14^ 

Xbitc) 
I 

eo 



dm 
dl 



dz \T<Tcut^t>' 



X : coordinate along the fiber target. 
y : coordinate orthogonal to the beam and fiber. 
z : coordinate along the COSY-beam axis 
tc : time during COSY's machine cycle 
number of COSY cycles measured 

thickness of the fiber in the vertical and longitudinal di- 
rections 

cross section of the fiber 
Avogadros number 
density of polypropylene 
molar mass of CH2 

horizontal and vertical beam width (standard deviation) 

horizontal beam position (centroid) 

beam current. 

elementary charge. 

beam momentum. 

mass per unit length of the fiber. 

restricted energy loss [7] of protons in polypropylene. 

nn : number density of hydrogen atoms in polypropylene. 

nH,o • number density of hydrogen atoms in undamaged 
polypropylene nH,o = SAT^^icHaZ-^CHa- 



density of hydrogen atoms normalized to the value of un- 
damaged polypropylene. 



QN = nn/nn,o 

D{x,t) : dose distribution along the fiber target as a function of 



linio. 



Table 1 

Notation and variables used throughout this paper. 



In section 2 the EDDA experiment is described, the formalism used to analyze 
the data in order to extract the hydrogen density as a function of dose is the 
topic of section 3. The results are presented and discussed in section 4. 
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Fig. 1. The EDDA detector (not to scale): target: fiber (CH2 or C); B: scintillator 
bars; R: scintillator semi-rings; F: semi-rings made of scintillating fibers. 



2 Experiment 



The EDDA experiment [1,8,9] has been designed to measure differential cross 
sections, analyzing powers and spin correlation parameters of elastic proton- 
proton scattering over a wide energy (0.5-2.5 GeV) and angular (30°- 90° in 
the cm.) range. A schematic view of the EDDA detector is shown in Fig. 1. 
It is a cylindrical scintillator hodoscope optimized to efficiently detect two 
protons from elastic pp scattering in coincidence. Elastic events are identified 
by coplanarity with the beam and the kinematic correlation of the scattering 
angles tan^iab,i • tan^^iab,2 = Ic.m.i where 7c.m. is the Lorentz-parameter of the 
cm. system in the lab. 

For this study scattering data acquired at fixed energies of 3.0 and 1.455 GeV/c 
with CH2 and carbon targets, detected primarily for the purpose of detector 
calibration, have been used. The fibers are suspended to a 30 mm wide fork 
and coated with 20 /xg/cm^ aluminum to avoid charge buildup on the target. 
The targets are mounted on a linear actuator and can be moved in and out of 
the beam. The experiment is conducted during many COSY machine cycles: 
after beam injection into COSY, the beam is accelerated to the desired beam 
momentum, then the fiber target is moved in and data is taken. Since beam- 
lifetimes are of the order of a few seconds the target is moved out after 5 to 
10 seconds, and COSY is reset to be prepared for the next cycle. 
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3 Analysis 



For the analysis data has to be sorted as a function of the time tc an event 
occurred within a cycle, i.e. tc is the time elapsed since moving in the target. 
Due to beam-heating by the target, the width and eventually the position 
of the COSY beam, and thus the dose distribution, change with tc, however, 
they are stable from one cycle to another. The target was first irradiated for 
about 5 h at a beam momentum of 3.0 GeV/c and then for another 5 h at 
1.455 GeV/c. The complete data set was subdivided in small sets comprising 
data of typically 70 machine cycles, corresponding to roughly 15 minutes of 
data taking. 



3. 1 Overview 



The dose per time interval dtc acquired over A^^^ machine cycles is a function 
of the horizontal position x along the fiber and the time tc during COSY's 
machine cycle. It is given by the product of the average energy deposit /S.E 
per particle and the number of particles hitting the target per unit length of 
the fiber divided by its mass per unit length: 



AE ANp/Atc beam profile 



d^("''^)= d7 



I(tc,c) dy^ 




(p) dz {J2 ^7-^ ) . (3) 



Of these quantities the beam current / and the vertical beam width ay are hard 
to measure precisely. However, the luminosity - with respect to hydrogen - 
along the target contains the same parameters 

AN p/ Ate beam profile 

target area density 



<1Lh, , ^I(tc,c)dy ■ — p 1 _(fz^j)! 

-^{^^'tc) = 2^ 7^= — nuflgN{x,tc,c)dz ^= — e _ 

ax ^^leo^/ In ay yl'nax 



In the experiment we access the luminosity by detecting elastic scattering 
events A^pp for A^^^ machine cycles as a function of tc- Due to the finite reso- 
lution in determining the x-coordinate of the scattering vertex only the total 
luminosity, e.g. integrated over x, is available: 

°^ AT N if \ 

LnitcNc)^ / dx-^{x,tc) = — • 

ax CTpp [1 - T{tc)\ e 
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Here, Upp is the differential cross section integrated over the accepted soUd 
angle, r the deadtime of the data acquisition system and e the detection 
efficiency, i.e. mainly a correction for losses of elastic events due to secondary 
reactions in the detector. 

Note, that the target density itself changes with the dose applied and will be 
both position and time dependent, so that the integration of Eq. (4) cannot 
be carried out directly. However, we can adopt an iterative approach and set 
qn to unity as a first approximation, which will overestimate the luminosity. 
The true luminosity will be smaller by the weighted average over x of the 
normalized hydrogen density, introduced via 



^/(tc, c) / dxgN{x,te,c) —= — e 

c=l _ V ^TT Cj; 

QN{tc) = — , (6) 

c=l 

a correction, we will have to determine using the relation Qn{D) to be estab- 
lished in Sec. 4, where in turn D is a function of x, tc and the individual cycle. 
Again in a first approximation pjv will be set to one. and recalculated in the 
following iterations. Now we may write Eq. (5) by inserting Eq. (4) as 

c=l ^0 V^TrCTy 



and thus 



2TTayeo n}i,oQN{tc) Na QcHi QN{tc) 



(8) 



Note, that this derivation holds, even when the vertical target position i/t is 
not centered at the vertical location of the beam yj,, since the additional factor 
exp{—[yt — yfjY/2ay) in Eqs. (3) and (4) cancels. 

Insertion in Eq. (3) finally replaces / and ay by experimentally accessible 
quantities 

k{p) 



dD^ , dE 
d^("''^)= d7 



This must be integrated over tc to obtain the dose profile D{x) along the beam 
acquired during Nc machine cycles. 
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Fig. 2. Distribution of tlie kinematic deficit a for data taken witli a CH2(sofid) and 
C-targets (open symbols, normalized to the CH2-data). The upper limit in a used 
to select elastic scattering events (cut) is indicated as well as the range (norm) used 
to determine the relative luminosity with respect to scattering of carbon for the two 
targets. 



The restricted enere;v loss 4^ [71 is the enersry loss due to ionization, 

corrected for the energy carried away by 6 electrons with energies T > T^ut 
escaping from the target. By requiring the practical range of electrons [10] 
to be half of the CH2-target thickness 0.2 mg/cm^) and taking the thin 
aluminum-coating (Ri 20 ^g/cvo?) into account we deduce a value of 9 keV for 
Tcuf Using ionization-loss parameters of [11] we obtain numerical values for 
the constant k{p) of 190.9 (155.0) Gymmmb at 1.455 (3.0) GeV/c. 



To summarize, we need to deduce the integrated luminosity Lh-, the beam 
position and width from the scattering data in order to obtain the dose profile 
along the fiber. As pointed out earlier, Eqs. (6) and (9) can only be solved 
iteratively. One starts by assuming ^(tc) being unity, and then arrives at a 
dose which will be underestimated. Prom the measured change in hydrogen 
density as a function of the dose (cf. Eq. (12)) ^(tc) is recalculated and a bet- 
ter estimate of the dose is found. After at most two iteration self-consistency 
is reached. 



3.2 Luminosity 



For the luminosity determination, using Eq. (5), we need the total number 
of elastic scattering events. These are identified by looking at the so-called 
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"kinematic deficit" , the angle a of the deviation from a perfect back-to-back 
correlation of the two detected prongs in the cm., when transformed assuming 
an angle-energy correlation as for pp elastic scattering. As evinced in Fig. 2, 
pp elastic scattering events stand out at very small a with very little back- 
ground from proton-carbon scattering. We use a cut at 10.2°(7.7°) at 1.455 
(3.0) GeV/c for accepted events. The geometric acceptance of the EDDA- 
detector is taken into account by using only events with polar angles between 
36° (38°) and 90° in the cm. for the forward proton. In addition, events with 
an azimuthal angle within ±15° around 90° or 270° are excluded, where the 
angle measurement of ^c.m. has reduced accuracy owing to the inhomoge- 
neous light-collection close to the attached light-guides. The contribution from 
proton-carbon scattering is removed by subtracting the yield from scattering 
of a carbon target using the same cuts and normalized to the same number of 
counts in the tails of the a-distribution as indicated in Fig. 2. 

The pp-elastic scattering cross section, cf. Eq. (5), is obtained by numerical 
integration of the differential cross section as given by the solution FAOO [12] 
of the VPI/GW phase shift analysis. 

Further corrections are applied to account for losses due to secondary reac- 
tions of the ejectiles in the detector (5%) and for the deadtime of the data- 
acquisition event-builder, deduced from the ratio of triggers read out by the 
DAQ-system to the total number of triggers, which is accomplished by 20 MHz 
scalers, read out every 2.5 ms. Details are given in [1,5]. The deadtime ranges 
from 90% at the beginning to 10% at the end of the cycle. The resulting 
luminosity is shown at the top of Fig. 3. 



3.3 Beam Parameters 



The beam position and width can be inferred from pp elastic scattering data as 
well, by exploiting the coplanarity of the two protons with the beam. The beam 
position and width as well as the azimuthal angular resolution are obtained 
from a nonlinear x^-^t of data with different values of the azimuthal angle 
4>. The procedure is described in detail elsewhere [6], the results for the two 
beam-momenta are displayed in Fig. 3. 

The increase of the beam-width or emittance due to the heating of the stored 
beam is clearly visible. At 1.455 GeV/c the movement of the beam at the end 
of the cycle is due to the beginning deceleration, whereas the jump in position 
at mid-cycle of the 3.0 GeV/c data is caused by exciting horizontal steerer 
used to cross-check the vertex-reconstruction. 
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Fig. 3. Luminosity (top), horizontal beam position (center) and width (bottom) as 
a function of time during the COSY-cycle (here tc=0 has been set to the time the 
target is moved in) for beam momenta of 1.455 and 3.3 GeV/c. 

3.4 Hydrogen Density 



The relative change of the hydrogen density in the target is proportional to 
the ratio of A^pp, i.e. the shaded area in Fig. 2, and the number of proton- 
carbon scattering events A^pc- For the latter we use the number of events in 
the range 11° < a < 25° of Fig. 2, mainly populated by quasi-free proton- 
proton scattering. To this end the normalized hydrogen density is given by 



Nr 



pC 



(10) 



with a momentum-dependent proportionality constant /(p) which needs to be 
determined from the data. 
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4 Results and Discussion 



The data has been subdivided into 43 subsets, 21 at 3.0 GeV/c and 22 at 
1.455 GeV/c, and a value for Qn/ f{p) has been deduced viz Eq. (10). This 
needs to be related to the dose acquired prior and during data taking of this 
subset. Since ^tv is averaged over the horizontal beam profile, we relate it to 
the corresponding average dose 



/ dte A^pp(tc) / dxD{x,Q -j= 
D^i , (11) 



where we imphcitly assume, that an appropriate summation over the num- 
ber of cycles has been done. For simplicity, we will denote this dose by D 
subsequently. 

For the functional dependence of ()n{D) we tried a simple ansatz of a linear or 
exponential decrease, which turned out to be insufficient to describe the data. 
We obtained a successful fit with the following expression: 

77,TT 

QNiD) = -^{D) = aexp(-AiL>) + (1 - a) exp(-A2£>), (12) 



using a total of five fit parameters a, Ai,2, /(1. 455 GeV/c), and /(3.0 GeV/c). 
The result is shown in Fig. 4. 

When iterating Eqs. (6) and (9), the average hydrogen density for a sample is 
given by 

D(,+AD 
Do 



where Dq and AD are the dose applied prior and during the measurement of 
the sample. After two iterations self-consistency is reached and we obtain the 
fit values 

Ai = (4.4 ± 1.3 ) -10-8 Gy-i 

A2 = (5.67 ± 0.28) -10-1° Gy-^ (14) 
a = (5.8 ± 0.7 ) % 
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Fig. 4. Normalized hydrogen density as a function of applied dose. The result of the 
fit is shown as the solid line. The dashed line shows the contribution of the slow 
component only. Note, that the first data point at 1.455 GeV/c corresponds to a 
lower dose than the last point at 3.0 GcV/c, since the beam is slightly displaced, 
such that a region of lower accumulated dose is probed. 



where errors are purely statistical. The systematic error is dominated by un- 
certainties in the area A of the target cross-section (10%) and the restricted 
energy loss (4%), where the exact value depends on the choice of Tcut- This 
leads to a 11% scale uncertainty of the dose and consequently of Ai and A2. 

The normalized hydrogen density qn of Eq. (10), determined from the experi- 
ment, is in fact the average over the length of the fiber. It receives contributions 
from regions with quite different accumulated dose, which are weighted with 
the beam intensity distribution and averaged over the machine cycle. It can 
be shown, however, that the correct functional dependence is obtained when 
Qn^D) is a linear function. The systematic eflPect caused by the deviation from 
hnearity in Eq. (12) was checked numerically and is always smaller than 0.25% 
in Qjsf, and can therefore be neglected. 

Assuming that the loss of hydrogen is entirely due to cross-linking, our results 
may be compared to the so-called G-value [4] found in the literature. It de- 
scribes the number of cross links SNx per 6E = 100 eV energy deposited in 
the polymer. Since per cross-link two hydrogen atoms are lost, the change of 
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Fig. 5. G-value for cross-linking as a function of the accumulated dose, 
hydrogen density for a sample of volume V and mass m is given by 

_d^_J_ 25Nx rn^^QcH^ 2g ^ Men, G_ 
dD ~ riH.o ' V ' SE nH,o ' SE Na ' SE' 



(15) 



Solving for G using Eq. (12) yields 

G = [aAi exp(-AiL') + A2(l - a) exp(-A2L')] (16) 



which is displayed in Fig. 5. It nicely agrees with the range of G-values 0.3 
to 1.1 [13] and 0.6 ± 0.1 [14] found in the literature. However, we find that 
the G-value is not a constant but depends on the irradiation-history of the 
target. The data of [14,15], obtained by a similar measurement with 14.5 MeV 
protons on 14 ^m thick polypropylene fibers but with much lower statistical 
precision, are in agreement with this finding. 



5 Summciry 



Thin CH2-fiber targets have been irradiated by 1.455 and 3.0 GeV/c pro- 
tons at an internal target station of the COSY accelerator. The applied dose 
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distribution was deduced from the elastic scattering rate measured continu- 
ously during irradiation. The hydrogen density was monitored by observing 
the relative yield of proton-proton elastic scattering with respect to proton- 
carbon inelastic reactions. The hydrogen loss as a function of dose can be 
described by a fast component responsible for at most 6% of the total loss, 
and a slow component corresponding to G-value of about 0.35. The loss rate, 
attributed to cross-linking of the chains within the polymer, is in-line with 
published values. Detailed knowledge of the observed strong dependence on 
the accumulated dose of a polypropylene-target will be important for experi- 
ments aiming at absolute cross-sections, where the hydrogen-density needs to 
be known accurately as a function of time as well as vertex position that may 
vary considerably during acceleration [5]. 
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